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Abstract 

In this paper we demonstrate that there exists a close relationship between 
quasi-exactly solvable quantum models and two special classes of classical dy- 
namical systems. One of these systems can be considered a natural generaliza- 
tion of the multi-particle Calogero-Moser model and the second one is a classical 
matrix model. 
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1 Introduction 



A quantum mechanical model is called quasi- exactly solvable (QES) if a finite number 
of energy levels and the corresponding wavefunctions of the model can be constructed 
explicitly. One possible way of studying the QES models and their solutions in the 
one-dimensional case is based on a reformulation of the spectral equations in terms 
of the wavefunction zeros. It can be shown that, in all the cases when the number of 
wavefunction zeros is finite (and this is just the case of QES models), the problem 
of reconstructing solutions of the stationary Schroedinger equation becomes purely 
algebraic and is reduced to the determination of their positions in the complex plane. 
These positions (which hereafter we shall denote by ^j) obey the system of algebraic 
equations 

M * 

E 73^ + ^(^0 = 0, i = l,...,M (1) 

where -F(^) is a rational function in which all the information of a QES model is 
contained. In turns out that equations of the type (|l|) are typical not only for QES 
models but appear in many branches of mathematical physics and this fact enables 
one to establish a deep relationship between QES models and many other seemingly 
unrelated models. For example, QES models turn out to be equivalent to com- 
pletely integrable Gaudin spin chains [l^ (for which system ([l|) plays the role 
of the Bethe ansatz equations for the coordinates of elementary spin excitations), 
to random matrix models [Q, |5[ (for which system (2) determines the distribution 
of eigenvalues of large random matrices), and also to purely classical models of 2- 
dimensional elctrostatics 0, ||, ^] respectively hydrostatics of point vortices . For 
both the latter models system (|l|) determines the equilibrium positions of pointwise 
classical objects (the charged Coulomb particles or resp. vortices) in an external 
(electrostatic or resp. hydrostatic) field. 

Up to now only the stationary solutions of QES Schroedinger equations have 

been considered from this point of view. The aim of the present paper is to consider 

time dependent solutions for QES models and to show how they can equivalently 

be described as dynamical equations for the motion of wavefunction zeros ^. This 

enables one to reveal three classical (complex) dynamical systems closely related to 

quantum QES models. One of these classical systems does not have an immediate 

physical interpretation, but the two remaining ones are very interesting from both 

the physical and mathematical point of view. The point is that one of these two sys- 

*We mean here just the solutions of evolution equations for standard QES models with time- 
independent potentials. Note that QES models with time-dependent potentials were discussed (from 
different point of view) in paper |7| 



2 



terns is a natural generalization of the famous classical Calogero-Moser multi-particle 
system and the second one is a classical matrix model from which the first system 
can be obtained by means of Olshanetsky and Perelomov's projection method ||l3|| . 



Stricly speaking, the idea of studying classical dynamical systems describing the 
motion of zeros of solutions of linear differential equations belongs to Calogero and 
is far from being new. Papers devoted to the investigation of such systems appear 
in the literature rather frequently. The goal of our paper is to apply Calogero and 
Olshanetsky-Perelomov methods to a concrete class of quantum QES models and to 
derive the associated classical models. The most interesting and quite unexpected 
result which we intend to present here is that the potentials of the resulting classical 
matrix models turn out to almost coincide with the potentials of the initial quantum 
QES ones (up to terms of order h) . This fact may hint to the existence of a certain 
non-standard "quantization procedure" relating the classical multi-particle systems 
of Calogero-Moser type to QES models of one-dimensional quantum mechanics. 

The paper is organized as follows: in Section 2 we remind the reader of the 
basic facts concerning the simplest QES model in the stationary case. In Section 3 
we derive the explicit form of solutions for this model in the time-dependent case 
and the corresponding system of evolution equations for the wavefunction zeros. 
In Section 4 we show thatsolutions of this system can be considered as "soliton 
like solutions" of a multi-particle classical system of the Calogero-Moser type. The 
matrix version of this system is discussed in Section 5. In Section 6 we discuss the 
limiting case when the number of algebraically calculable states in the QES model 
tends to infinity. The last section is devoted to a discussion of our results. 



The simplest QES model is the sextic anharmonic oscillator. Its potential has the 
following form 



where a and b are real parameters with a > 0, and M is an arbitrary non-negative 



2 The simplest QES model 







admits solutions of the form 



M 




(4) 



i=l 
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E = hbi^M + -j+haY,^f, (5) 

where the numbers ^j, i = 1,...,M (playing the role of the wavefuntion zeros) 
satisfy the system of numerical equations 

M * 

E F3^ = ^^^ + <'' ^ = l>--->^ (6) 



with the following additional condition 



M 



i=l 



Note that the form of system (^) coincides exactly with that of the famous Bethe 
Ansatz equations appearing in the theory of completely integrable Gaudin models 
pi IC]. Therefore we hereafter will call (pi) the Bethe Ansatz equations. 



It is not difficult to show that for any fixed M the Bethe Ansatz equations (P)-(|7|) 
have 



M 



+ 1 (8) 
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solutions. This means that for the sextic anharmonic oscillator we can find Ng^i 
wavefunctions ^i{x) and N^oi corresponding energy levels Ei by means of purely 
algebraic methods. 

Note also that model (^) can be considered a deformation of the simple harmonic 
oscillator with the potential 

Vo{x) = (9) 

The role of the deformation parameter is played by a. It is not difficult to see 
that after taking a = in formulas (0)-(|6|) they reduce to the well known formulas 
describing the solution of the simple harmonic oscillator in terms of wavefunction 
zeros Ijl], In this case the number M can be considered a free parameter because 
it is not any longer correllated with the form of the potential. 

3 Evolution equations 

Let us next consider the time-dependent Schroedinger equation 

8^(x,t) i h'^ ,] , 

^^4^= --^ + y(x)U(x,t). (10) 
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for potential (|2|) . It is obvious that any linear combination of the stationary solutions 



.n=l i=l 



n=l 



exp 



1 / ax^ 



ill) 



gives a certain dynamical solution of equation (p!o|). Remember that the pre-exponential 
factor in ^ is always a polynomial of degree M. Therefore also the pre-exponential 
factor in (^) is a certain polynomial and hence its zeros are functions of time. This 
enables one to write 



M 



^{x,t) = Cit)l[ix-^i{t))exp 



i=l 



1 / bx"^ ax^ 



(12) 



Formula (^) enables one to derive the evolution equations for the functions ^i{t) 
and C{t). For this it is convenient to rewrite equation (Il0|) in the form 



d 



dt 



V{x) = ifi— In ^{x, i) + — i ( ^ In ^{x, t)] + ^ In ^{x, t) 



dx 



(13) 



Substituting expressions (§) and (|^) into equation ( [l^ ) we obtain after some algebra 
the relation 



C{t) 



M M / 1 \ 

E iKt) -haxJ2m-f^b(M + -) + 
i=l i=i \ z/ 



M 



M 



E 



^m- E 



+ bm + aCfit)\ =0 (14) 



from which it immediately follows that the functions £,i{t) must satisfy the following 
system of first order differential equations 



M 

iUt)= E 



k=l,k^i 

supplemented by the condition 



bUt)-<fit), i = l,...,M 



(15) 



M 



E^^w = o- 



(16) 



i=l 



For the function C{t) we obtain on the other hand 



Cit) = exp 



-I a 



E / ^f{t)dt + b(^M + ^^ 



(17) 
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The system ([T5|)-(|l6|) is the dynamical extension of the stationary system 
Since the functions ^j(t) are assumed to be complex, we face a typical example of a 
complex dynamical system of dimension M. It is remarkable that this system can 
be rewritten in the following "potential" form 

iUt) = -Q^um), (18) 

where 

M , M M 

u{i) = -n Hm-ut)) + \Y.iKt) + ^^Y.it{t) (19) 

k=l,kj^i i=l i=l 

is playing the role of a complex potential. 



4 Complex multi-particle systems 



Despite the fact that system ( piq ) cannot be directly derived from the Lagrange 
principle, it is possible to relate it to a certain Lagrangian system in the following 
way. For this, consider the complex "Lagrangian" 



1 ^ 1 



with 



i=l 



(20) 



(21) 



It is not difficult to see that this Lagrangian (^) can be represented in the form 

M 



i=l 



where 



M 



n 



(22) 



(23) 



k=l,k^i 

The total time derivative in (|2^ ) can be omitted because it does not affect the form 
of the equations of motion. The form of these equations derived from the action 
principle reads then 



M 

E 

i=l 



d 
It 



din 



din 



0, n = l,...,M. 



From (24) it immediately follows that the solutions 

Ai(e,O = 0, i = l,...,M 



(24) 



(25) 
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are automatically solutions of the dynamical equations for the Lagrangian (^) . Note 
that system (|25| ) exactly coincides with equations (|l5|). To derive the form of the 



Lagrangian (^) it is sufficient to substitute expression (|19|) into (|20|). This gives 

, M M M ^2 

- E ^ + ^) + - E (26) 

We have obtained a complex Lagrangian system supplemented with the additional 
constraint (|l^). This system can obviously be regarded a deformation of a complex- 
ified version of the famous Calogero-Moser system with Lagrangian 

1 M M M *2 t2 M 

Note that the above construction is very similar to the one usually used in con- 
structing the soliton solutions of some classical dynamical equations. Therefore it 
is natural to interpret the solutions of equation (|2^) describing the motion of wave- 
function zeros in quantum QES model (|2|) as solitons of the classical multi-particle 
system with Lagrangian (|20|). The relation between the undeformed theories of 
the simple harmonic oscillator and the Calogero-Moser system is well known in the 



literature |11, 12 1. 



5 Classical matrix model of the sextic anharmonic os- 
cillator 

Let X be a M X M traceless complex matrix whose entries are considered as dy- 
namical variables of a certain complex dynamical system. The Lagrangian of this 
system can be chosen in the form 

C{X, X) = ^TrX^ - TyV{X) (28) 

with V{X) given by formula (Q). The corresponding dynamical equations then read 

X = -^Vin (29) 

It is known that any complex matrix with non-coinciding eigenvalues can be reduced 
to diagonal form by means of an appropriate similarity transformation 

X = SES-^ (30) 
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where we used the notation E = diag{^i(t), . . . ,S,Ai{t)} with ^i(t) + • • • + (,M{t) = 0. 
Taking the time derivative of ( |30|) we obtain 



X = SLS~^, (31) 

where 

L = E + [M,E] (32) 

and 

M = S-^S. (33) 
Differentiating (31) once more and using (p9[) and (pO) we get the equation 



L + [M,L] = -^V{E) (34) 
which together with (^2|) is equivalent to system (p9|). Note that (^) has the form 



of a deformed Lax equation. A relation between equations (|2^) and (^4^ can also 
be established in the opposite direction: Assume that we have two matrices M and 



L satisfying (pj4|) and (|3^). Let 5 be a solution of the linear evolution equation 

S = SM. (35) 



Then the function X in ( pO| ) is a solution of equation ([29|). It is easy to check that 
the matrices L and M with components 

_ ■ , , i/i 



satisfy the system (p2D, ( p4\) provided the functions ^i{t) are solutions of the deformed 
Calogero-Moser system ([2^). As mentioned above for this it is sufficient that (,i{t) 
be solutions of system (^5[)-(^) describing the motion of the wavefunction zeros in 
the quantum QES model with potential (|2[). 

From the above reasonings it follows that we essentially established a relationship 
between the quantum QES model with Hamiltonian 



2 2 

and the classical matrix model with Hamiltonian 



H = ^ + — fax^ + b) -ha(M + -)x^ (37) 
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6 The large M limit 



It is remarkable that the hamiltonians of both the quantum and classical models 
( |37| ) and ( |38| ) essentially coincide. The difference is of order ha. Let us demonstrate 
now that this difference also disappears if we take a proper large M limit of these 
models. 

The neccessity for taking the limit M ^ oo in formulas (^7|) and ( |38| ) arises from 
the following reasonings. For finite M it is not correct to speak of an equivalence 



of models ( pTf ) and (38). The point is that for finite M model ( |37|) describes the 
evolution of only a certain finite superposition of quantum states (because the model 
( |37| ) is quasi-exactly solvable). However, if M tends to infinity, then the number of 
stationary states in this superposition also tends to infinity and fills all the spectrum 



of the model. Only in this case we can say that models (37) and ( |38[) are equivalent. 
But how to proceed to the large M limit? It is clear that if we simply take M = oo 
in formulas (^^ and ( |38| ) then we obtain a meaningless (minus) infinity. In order to 
get a finite expression we must take into account that the parameters a and b may 
also depend on M. Choosing this dependence according to the conditions 

ab = g, — - haM = — (39) 

in which g and oo are fixed (positive) numbers, we find that in the large M limit 
the parameter a must behave as a ~ M~^l^. This means that in the limit M — > 00 
the terms containing a? and a (i.e. the sextic term and also the harmonic term 
responsible for the difference between the models) disappear and we obtain the two 
models 

F=^ + i^ + ,.^ (40) 

and 

W = Tr(^ + i^+,X^) (41) 

with exactly coinciding classical and quantum hamiltonians. Note however that 
model (^) is a one-particle quantum model while model (^) is a classical model of 
traceless infinite matrices. The models are equivalent in the sense that the soliton- 
like solutions in the classical model (pl| ) describe the evolution of the wavefunctions 
in the quantum model (|40| ) and vice verso. 



7 Discussion 



When speaking of a relationship between quantum and classical mechanics one usu- 
ally thinks of a pair of quantum and classical models related to each other by a 
certain quantization - dequantization procedure. The hamiltonians of these models. 
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considered as functions of coordinates and momenta, formally coincide up to terms 
of order h. However, the mathematical meaning of the coordinates and momenta 
is essentially different in the quantum and classical case. In the quantum case they 
are operators in Hilbert space (infinite matrices) while in the classical case - they 
are simply numbers. Correspondingly, the number of quantities (degrees of free- 
dom) neccessary to fix uniquely the state of a dynamical system is infinitely larger 
in the quantum case compared to the classical one. In this sense the transition to 
the classical limit is equivalent to freezing infinitely many degrees of freedom of the 
quantum system. 

There are however several examples where the correspondence between the quan- 
tum and classical model is not approximate but exact and is not related to any 
limiting procedure and any lose of degrees of freedom. The idea underlying such 
examples is based on a proper parametrization of the wavefunction by an infinite 
number of parameters depending on time and considered as canonically conjugated 
dynamical variables. If such a parametrization is found then we can associate with 
a given quantum model a certain classical one which in this case should neccessarily 
be infinite-dimensional. It is quite clear that a priori there are no reasons for any 
relation between the forms of the corresponding quantum and classical hamiltoni- 
ans. In general, they may be of absolutely different nature. Consider a simple but 
most instructive example. Let us take the evolution equation for a certain quantum 
model with hamiltonian H 

ihdt^ = i^^- (42) 

Since the wavefunction is complex, "if = Q + iP , we can rewrite equation (^) in 
real form 

HdtQ = HP, hdtP = -HQ, (43) 
which, after introducing the functional 

n{Q, P) = ^ {Q, HQ) + ^ (P, HP) , (44) 

can in turn be rewritten in the form of the classical Hamilton-equations: 

a,«=^W§Z) 3.P=J-m^. (45) 

oP OQ 

We see that irrespective of the form of the initial quantum model, the resulting 
classical one describes an infinite-dimensional coupled harmonic oscillator. 

In this paper we have found examples for an exact relationship between quantum 
and classical models. The main distinguished feature of these examples is that the 
potentials of the initial quantum and the resulting classical models exactly coincide. 
In this case the construction of the classical counterpart of a given quantum model (in 
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our case the quartic or harmonic oscillator) is extremely simple. One should simply 
replace the operators of coordinate and momentum entering into the quantum model 
by infinite traceless complex matrices and, after this, take the trace of this matrix 
hamiltonian. What we obtain will be just the hamiltonian of a classical model whose 
solutions contain the complete information of the dynamics of the wavefunctions in 
the initial quantum model. The origin of this coincidence is not clear to us at the 
moment but it is quite obvious that it cannot be accidental. It would be tempting 
to conjecture that any one-dimensional quantum model with hamiltonian 

H = ^ + V{x) (46) 

is somehow equivalent to the classical infinite and traceless complex matrix model 
with the same hamiltonian 



n = TT\— + ViX) (47) 




We have an idea how to check this conjecture at least for models with polynomial 
anharmonicity and hope to publish the results in the near future. 
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